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Abstract 

We study a well-known problem concerning a random variable Z uniformly distributed 
between two independent random variables. Two different extensions, conditionally directed 
power distribution and conditionally undirected power distribution, have been introduced for 
this problem. For the second method, two-sided power random variables have been defined. 

1 Introduction 

Van Assche (1987) introduced the notion of a random variable Z uniformly distributed between 
two independent random variables X 1 and X 2 , which arose in studying the distribution of products 
of random 2x2 matrices for stochastic search of global maxima. By letting Xi and X 2 to have 
identical distributions, he derived that: (i) for X\ and X 2 on [— 1, 1], Z is uniform on [—1, 1] if and 
only if Xi and X 2 have an Arcsin distribution; and (ii) Z possesses the same distribution as X\ 
and X 2 if and only if X 1 and X 2 are degenerated or have a Cauchy distribution. Soltani and Homei 
(2009) following Johnson and Kotz (1990) extended Van Assche's results. They put Xi, ■ ■ • , X n 
to be independent, and considered 

S n = RiX\ + R 2 X 2 + • • • + R n -iX n -i + R n X n , n > 2 , 

where random proportions are R4 = U^) — Lfy-i), % — 1, n — 1 and R n — 1 — Y17=i R%, 
U(i), U( n -i) are order statistics from a uniform distribution on [0, 1], and U(p) = 0. These ran- 
dom proportions are uniformly distributed over the unit simplex. They employed Stieltjes trans- 
form and derived that: (i) S n possesses the same distribution as Xi,...,X n if and only if Xi,...,X n 
are degenerated or have a Cauchy distribution; and (ii) Van Assche's (1987) result for Arcsin holds 
for Z only. 

In this paper, we introduce two families of distributions, suggested by an anonymous referee 
of the article, to whom the author expresses his deepest gratitude. We say that Z\ is a random 
variable between two independent random variables with power distribution, if the conditionally 
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distribution of Z\ given at X 1 — x\ , X 2 — x 2 is 

1 z > max(xi, x 2 ), 



F Zl \xi,X2(z) 



(*=2±y Xl <z< x 2 , 

z < min(xi, x 2 ). 

The distribution F^Uj^ will be said to follow a conditionally directed power distribution, when 
n is an integer. For n = 1, the distribution given by (1.1) simplifies to the distribution Z that was 
introduced by Van Assche (1987). For n = 2, we call Z x directed triangular random variable. 
For further generalizing Van Assche results, we introduce a seemingly more natural condition- 
ally power distribution. We call Z 2 two-sided power (TSP) random variable, if the conditionally 
distribution of Z 2 given at X\ — x± , X 2 — x 2 is 

1 z > y 2 , 

Fz 2 \ xl , X2 (z) = { (J^)™ yi<z<y 2 , (1-2) 
z< Vl . 

The distribution F Z2 \ xltX2 will be said to follow a conditionally undirected power distribution, when 
yi = min(xi, x 2 ), y 2 = max(xi,x 2 ) and n is an integer. For n = 2, we call Z 2 undirected 
triangular random variable. 

Again, for n = 1, the distribution given by (1.1) simplifies to the distribution Z that was introduced 
by Van Assche (1987). The main aim of this article is providing a couple of generalizations to the 
results of Van Assche (1987) for some other values of n (other than n = 1). This article is organized 
as follows. We introduce preliminaries and previous works in section 2. In section 3, we give some 
characterizations for distribution Z x given in (1.1), when n = 2. In section 4, we find distribution 
of Z 2 given in (1.2) by direct method, and give some examples of such distributions. 



2 Preliminaries and previous works 

In this section, we first review some results of Van Assche (1987) and then modify them a little bit 
to fit in our framework, to be introduced in the forthcoming sections. 
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Using the Heaviside function (U(x) = 0, x < 0, =1, x > 0) we conclude that for any given 
distinct values x\ and x 2 , the conditional distribution F Zl \ Xl:X2 (z) in (1.1) is 

Fz llxi , x M = (^^) n U(z - Xl )-j2 ( U ) (^^TU(z - x 2 ). (2.1) 
x 2 - Xi \l J x 2 — Xi 

1=1 x ' 

Lemma 2.1. For distinct reals x\, x 2 , z and integer n, we have 

-1 (-1)™ d n ~ x . 1 1 1 
1 — - — ( . ) = . 

(z — xi)(x 2 — xi) n (n - 1)\ dx^' 1 z - x 2 (xi - x 2 ) (xi - z)(x 2 - z) n 



Proof. It easily follows from the Leibniz formula. □ 
Another tool for proving our main theorem is the following formula taken from the Schwartz 
distribution theory, namely, 

^x)A^(dx)= [ —J-j —<p(x)A(dx), (2.2) 

where A is a distribution function and A^ is the n-th distributional derivative of A. 
The conditional distribution F Zl \ Xl ,x 2 ( z ) given by (1.1) leads us to a linear functional on complex- 
valued functions / : R — > C, defined on the set of real numbers M: 



Jn—t 



rz ^ AJ) ( X2 - Xl )n ^ (n - i)\(x 2 - Xi y dz^ J[X2) - 
It easily follows that 

F Zl \ Xl ,x 2 (af + bg) = aF Zl 

\Xl,X2 (/) + bF Zl (<?), (2-3) 

for any complex- valued functions /, g and complex constants a, b. We note that F Zl \ xliX2 (z) = 
F Zl \ XUX2 (f z ), whenever f z (x) = (z - x) n U{z - x) and 

f z (x±) n ^ i d n ~ i 

F Zl \ Xl , X2 {fz) = (a?2 Z _ xi)n ~ g {n-i^-x^dz^^' 

Also we note that U(z - x) = { -^£^fz{x). Thus 

2 

P(^i<*) = / U{z-x)dF Zl {x)= [ F Zl \ xuxa (z)f[F x .(dxi), 
Jr Jr 2 i=1 



can be viewed as: 



(-l) n d 
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R ( n V- dx- 



■n 



f z (x)dF Zl (x) 



F Zl \x u x 2 (fz)Y[Fx i (dx i ). 



(2.4) 



i=l 



Therefore by using (2.3) along with (2.4) and a standard argument in the integration theory, we 
obtain that 



for any infinitely differentiable functions / for which the corresponding integrals are finite. Now 
(2.5) together with (2.2) lead us to 



for the above mentioned functions /, where F^' is the (ra)-th distributional derivative of the dis- 
tribution of Z\. 

Let us denote the Stieltjes transform of a distribution H by 



for every z in the set of complex numbers C which does not belong to the support of H, i.e., 
z E C n (suppH) c . For more on the Stieltjes transform, see Zayed (1996). 

The following lemma indicates how the Stieltjes transform of Z\, and Xi, X 2 are related. 

Lemma 2.2. Let Z\ be a random variables that satisfies (1.1). Suppose that the random 
variables X 1 and X 2 are independent and continuous with distribution functions F Xl and F X2 , 
respectively. Then 




(2.5) 




(2.6) 





2 



i=l 



Proof. It follows from (2.6) that 
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i=l 



and 
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for g z (x) = j^. Now, it follows that 

Fz\\xi,x2^9z) 7 rr / 7 TTT 

I .To — 7?1 I Z ' [ n, — 1. 1 1 Tn — 



d" • 1 



(#2 — ^i) n (n — i)!(x2 — xi) 1 dz n 1 z — x 2 

and by using Lemma 2.1, we have 

(_l)n 

Fzi\xi,x2 \9z) 7 T7 rr- 

(z — — cc 2 ) 

Therefore, 
and 

i 2 

-S^(F Zl ,z) = -S(F Xl ,z)S^- 1 \F X2 ,z), z e Cf)(suppF Xi ) c . (2.7) 

This finishes the proof. □ 
Note that Van Assche's lemma is the case of n = 1: 

-S'(F Zl ,z)=S(F Xl ,z)S(F X2 ,z). 

We also note that the Stieltjes transform of Cauchy distribution, i.e., S(F, z) = satisfies (2.7). 



3 Directed triangular random variable 

Let us now review Van Assche's result for directed triangular random variables. 

Theorem 3.1. If Xi and X 2 are independent random variables with a common distribution F x , 
then the characterizations of Z 1 for n = 1 and n = 2 are identical. 

Proof. We note that Xi and X 2 have a common distribution function F x . By using Lemma 
2.2 for n = 2, we have 

~S"(F Zl ,z)=S(F x ,z)S'(F x ,z), 

and so 

-S"(F Zl ,z) = ^-S 2 (F x ,z), 
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and 

-S'(F Zl ,z)=S 2 (F X: z). (3.1) 

We note that the Stieltjes transform tends to zero, when z is sufficiently large. In that case the con- 
stant in the differential equation will be zero. The equation (3.1) is exactly the equation obtained 
by Van Assche (1987) when X 1 and X 2 have a common distribution; so his results hold in our 
framework as well. □ 
This clever proof is due to the anonymous referee. Now, we apply Lemma 2.2 for some character- 
izations, when X 1 and X 2 are not identically distributed. 

Theorem 3.2. Let Xi and X 2 be independent random variables and Z\ be a directed triangular 
random variable satisfying (1.1). For n = 2, we have, 

(a) if Xi has uniform distribution on [—1,1], then Z x has semicircle distribution on [—1, 1] if 
and only if X 2 has Arcsin distribution on [—1,1]; 

(b) if X 1 has uniform distribution on [—1,1], then Z\ has power semicircle distribution if and 
only if X 2 has power semicircle distribution, i.e., 

m = 3 -^^, -1<z<1; 

(c) if Xi has Beta(l, 1) distribution on [0,1], then Z 1 has Beta(|, |) distribution if and only if 
X 2 has Beta(|, |) distribution; 

(d) if X 1 has uniform distribution on [0, 1], then Z\ has Beta(2, 2) distribution if and only if X 2 
has Beta(2, 2) distribution. 

Proof, (a) For the "if part we note that the random variable X\ has uniform distribution and 
X 2 has arcsin distribution on [-1,1]; then 

S{F x „z) = )fi*\z + \\-\*\z-\\), 

and 

S(F X2 ,z) 1 



Vz^T 

From Lemma 2.2 and substituting the corresponding Stieltjes transforms of distributions, we get 
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S (F Zl ,z) = 



V-l)§ 



The solution S(F Zl , z) is 

S(F Zl ,z) = 2(z-vG^l), 
which is the Stieltjes transform of the semicircle distribution on [—1, 1]. 

For the "only if part we assume that the random variable Z x has semicircle distribution. Then it 
follows from lemma 2.2 that 

S(F X2 ,z ' 



1-z 2 (^ 2 -l)i 
The proof is completed. 

(b) By an argument similar to that given in (a) and solving the following differential equations, 

S"(F z ,z) = 2 + z 2 )(\n\z + l\ -\n\z-l\)), (for the "if part), and 

(z z — 1) 2 4 

1 o/r. \ 3 2z+( y l-z 2 )(\n\z + l\-\n\z-l\) 

S(F X2 ,z) = i i ! 1Z 5 (for the "only if part), 



1-z 2 v ' ' 4 (l-z 2 ) 
the proof can be completed, 
(c) By Lemma (2.2), we have 



~S"(F z ,z) = — — L^_=L==, (for the "if part), and 

2 z{z-l) y/ z (z-l) 

' S(F X2 , z) = ~ ) , (for the "only if part) . 



- 1)^(3 - 1) 

The proof can be completed by solving the above differential equations, 
(d) By Lemma (2.2), we have 

S"(F Zl ,z) = — 7\( 6 ( z2 ~ z )( ln \ z \ ~ ln \ z ~ !|) - 6^ + 3), (for the "if part), and 

z(z — 1) 

S(Fx 2 ,z) = 6(z - z 2 ){ln \z\ - \n\z- 1|) + 6* - 3, (for the "only if part). 
Solving the differential equations, can complete the proof. □ 
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4 TSP random variables 



In section 3, we used a powerful method, based on the use of Stieltjes transforms, to obtain the 
distribution of Z x given in (1.1). It seems that one can not use that method to find the distribution 
of Z 2 given in (1.2). So we employ a direct method to find the distribution of Z 2 . 

Let us follow Lemma 4.1 to find a simple method to get the distribution of Z 2 . The work of 
Soltani and Homei (2009b) leads us to the following lemma. 

Lemma 4.1. Suppose W has a power distribution with parameter n,n> 1, n is an integer, and 
let Y\ = Min(Xi, X 2 ), Y 2 = Max(Xi, X 2 ), where X\ and X 2 are independent random variables. 
Let 

X = Y 1 + W(Y 2 -Y 1 ). (4.1) 

Then 

(a) X is a TSP random variable. 

(b) X can be equivalently defined by 

x = l -{x 1 + x 2 ) + (w- l -)\x, - X 2 \. 

Proof, (a) 

F x \ XuX2 (z) = P(Y 1 + W(Y 2 -Y 1 )<z\X 1 = x 1 ,X 2 = x 2 ) 
= P(yi + W(y 2 - yi )<z) 

2/2-2/1 

(b) The proof can be completed by substuting Min(X 1 , X 2 ) and Max(X l5 X 2 ) with Yi and Y 2 
in (4.1). □ 

4.1 Moments of TSP random variables 

The following theorem provides equivalent conditions for f/ k = EZ 2 . 

Theorem 4.1.1. Suppose that Z 2 is a TSP random variable satisfing (1.2). If X 1 and X 2 are 
random variables and E\Xi\ k < oo, % — 1, 2, for all integers k, then 



^\ T?7k _ „ r (fc+ 1 ) T(k-i+n) T?(\ri\rk-i\. 

{&) - n r(fc+n+1) 2^ i=0 r(k-i+l)^y Y l Y 2 J> 

(b) EZ\ = Eto 0(\) k ~ l E{W - + X 2 ) fc - l |X! - X 2 \ i ; 

(c) £Z 2 fe = Eto O^EiYi^Yi ~ YiY)- 
Proof, (a) By using Lemma 4.1, we obtain that 

EZ\ = E(J2 ( k ^j (1 - W) i YtW k - i Yf- i ) 

= ( k ^j E(W k -%l - W) i )E(YjY 2 k - i ) 
i=o 

_ n r(fc + i) y ^-. + n) 

- n r(A: + n+l)^r(A:-2 + l) 1 2 " 



(b) This can be easily proved by Lemma 4.1(b). 

(c) It straightforwardly follows from (4.1). □ 
Let us consider expectation and variance of Z 2 . First, we suppose that EYi = EY 2 = n 2 , 

VarYi = af, VarY" 2 = a 2 , and Cov(Y"i, Y" 2 ) = a 12 . Then 

„ „ = /ii + nn 2 
2 n+1 ' 

and also, if EX ± = EX 2 = 0, then 

71 

E(Z 2 ) = EYi + -(EY 2 - EYi). 

TX | X 

By Y 1 + Y~ 2 = X 1 + X 2 , we have 

T? 1 — 77/ 

E(Z 2 ) = E(Y 1 ) + ——(-2^) = T —EY 1 . (4.2) 

It can easily follow from (4.2) that the Arcsin result of Van Assche (1987) is only true for n — 1, 
and also, one can see that Theorem (3.2) in section 3 does not hold for the above Z 2 . 
About the variance, we have 

n({ii - /i 2 ) 2 + n(n + \fa\ + 2(n + l)(a 2 + na 12 ) 



VarZ 2 = 



(n+l) 2 (n + 2) 



10 



Following the computation of expectation and variance, we evaluate them for some well-known 
distributions. If Xi and X 2 have standard normal distributions, then from Theorem 4.1.19b) and 
the fact that X 1 — X 2 and Xi + X 2 are independent, it follows that their first, second and third 
order moments are equal, respectively, to 

V7r n + 1 

r,ry2 n 2 + n + 2 
EZz = — — -, and 



(n + l)(n + 2) 
i on~ -|- iz/; 2 

7 2 



P73 1 5n J + 12n^ + 13n - 30 



2^ (n + 3)(n + 2)(w + l) ' 
Also, in case X x and X 2 have uniform distributions, Theorem 4.1.1(b) implies that, 

k 

i ( /•• - I ) 

EZ k = n- 

i i) + /,• + I : 

i=0 



r(fc + l) >^ r(A:-i + n) 2 
" ' ~ n T(n + k + 1) ^ r(A; - i + 1) (A; + 2)(i + 1) ' 



2n+l 
£^2 = and 
3(n + 1) 

TA , 1 n 3 + 3n 2 + 6n + 2 

VariZ2) = 18 (n + l) 2 (n + 2) " 
Since some distributions do not have any moments, Theorem 4.1.1 is not applicable for investigat- 
ing Van Assche's results for them, whence, we prove the following theorem: 

Theorem 4.1.2. Suppose that Z 2 is a TSP random variable satisfying (4.1). Then 

(a) Z 2 is location invariant; 

(b) if X 1 and X 2 have symmetric distribution around /i, then Z 2 has symmetric distribution around 
/j,, only when n — 1 . 

Proof. 

(a) Is immediate. 

(b) We can assume without loss of generality that fj, — 0. If Z 2 has a symmetric distribution 
around zero, then 

Y 1 + W(Y 2 - Y 1 ) I -[Y, + W(Y 2 - Y ± )}. 



We note that 



Y\ + W{Y 2 - Y ± ) = [-Y, + W(-Y 2 - (-Y,))}. 
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Since, -Mm{X 1 ,X 2 ) = Max(-X l5 -X 2 ), X 1 = -X 1 and X 2 = -X 2 , we have 



Y 1 + W(Y 2 - Fx) = Y 2 + W{Y l - Y 2 ). 



(4.3) 



By equating the conditional distributions given at X x = x\ and X 2 = x 2 in (4.2), we conclude that 



It can also easily follow from Theorem (4.1.1) that the Cauchy result of Van Assche (1987) is true 
only for n = 1 . 

4.2 Distributions of TSP random variables 

In this subsection, we investigate computing distributions by the direct method. We will give two 
examples of derivation based on (4.1). This method may be complicated in some cases, but we 
have chosen some easy to follow examples. 

Example 4.2.1. Let Xi,X 2 and W be independent random variables such that X 1 and X 2 are 
uniformly distributed over [0,1], and W has a power function distribution with parameter n. We 
find the value fz 2 (z; n) by means of fz 2 \w(z\w); therefore 



By using the distribution of W, the density function f z , 2 (z; n) can be expressed in terms of the 
Gauss hypergeometric function F(a, b, c; z), which is a well-known special function. Indeed ac- 
cording to Euler's formula, the Gauss hypergeometric function assumes the integral representation 



where a, b, c are parameters subject to — oo < a < +oo, c > b > 0, whenever they are real, and 
z is the variable (see Zayed 1996). By using Euler's formula, the density function of Z 2 can be 
expressed as follows: 



n 



= 1. 



□ 




(4.4) 




fz 2 (z] n) 



2nz 



(1 - z n ~ l ) + 2(1 - z)z n F(l, n,n+l,z), < z < 1, 



(4.5) 



n-\ 
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where n > and n^l. When n — 1, similar calculations lead to the following distribution 



The probability density function fz 2 ( z ) was introduced by Johnson and Kotz (1990), for the 
first time, under the title "uniformly randomly modified tine". So f z . 2 (z;n) can be seen as an 
extension of the above mentioned distribution. We note that, from (4.1) and a simple Monte Carlo 
procedure using only simulated uniform variables, one is able to simulate the distribution (4.5). 

Example 4.3.1. Let X± and X 2 be independent random variables with Beta(l, 2) distribution. 
Then if W has Beta(3, 1) distribution, Z\ has Beta(2, 3) distribution. 

In the following theorem we compute the Stieltjes transform of Z 2 forn = 2. Let us remark 
that the complexity of the integral in the theorem indicates that for this case the direct method is 
preferred. 

Theorem 4.4.1 Let Z 2 be a undirected triangular random variable that satisfies (1 .2). Suppose 
that the random variables X\ and X 2 are independent and continuous with the distribution functions 
F Xl and F Xi , respectively. Then 



fz 2 (z) 



2(1 - z)log(l - z) - 2zlog(2), 0<z<l. 



-S'" (F z ,z)=S' (F Xl , z)S' (F X2 , z) + 2S(F Xl , F x , , z) 



where 



Proof. By using an argument similar to that given in Section 3, we can conclude that 




So, 




for# 2 (x) 



. From 



Fz 2 \xi,x2(9z) 



(z-xi) 2 (Z~X2) 2 
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and by using partial fractional rule, we have 

1 2 1 

Z ^' xM ~ ( Z - Xl )2( Z - X 2 f + ( X2 - Xl )2( Z - Xl )( Z - X2 y 

Therefore, 

2 



and 



~S"' (F Z2 ,z)=S' (F Xl , z)S' (F X2 , z) + 2S(F Xl , F x , 2 , z) . 



This finishes the proof. □ 
It is worth mentioning that the present method yields other extensions too; the following is 
such an example. 

Example 4.3.2. Suppose that Xi, X 2 , W are independent random variables. If X\ and X 2 have 
uniform distributions on [0, 1] and W has Beta(2, 2) distribution, then Z 2 has the same distribution 
as W. 

If the product moments of order statistics are known, those of W can be derived from that of 
Z 2 by using Theorem 4.1.1(a). Then the distribution of W is characterized by that of Z 2 . 

By an argument similar to the one given in Example 4.2.1, when W has a Beta distribution 
with parameters n and m, we find the distribution fz 2 { z \ n > m ) as 

V^2z 1 - I z (n - 1, m + ^21 - * J 2 (n, m - 1 , < z < 1, 

B(n,m) B{n,m) 

where I x (a, b) is incomplete Beta function: 

I x (a,b) = -^-r- I* t^il-tf-'dt, (a,b>0). 
B(a,b) J Q 

5 Conclusions 

We have described how (a) methods of Stieltjes transform, and (b) directed methods, could be used 
for obtaining the distributions, characterizations and properties of the random mixture of variables 
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defined in (1.1) and (1.2). Of course each one of the methods (a) or (b) has its own advantages 
and disadvantages, and none of them has a preference over the other. The TSP random variable 
when Xi and X 2 have uniform distributions, led us to a new family of distributions which can be 
regarded as some generalization of "uniformly randomly modified tine". The proposed model in 
the direct method can easily lead to distribution generalizations, though this is not possible for the 
first method, but here the characteristics can be easily computed. 
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